Internal state conversion in ultracold gases 
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We consider an ultracold gas of (non-condensed) bosons or fermions with two internal states, and 
study the effect of a gradient of the transition frequency between these states. When a 7r/2 RF 
pulse is applied to the sample, exchange effects during collisions transfer the atoms into internal 
states which depend on the direction of their velocity. This results, after a short time, in a spatial 
separation between the two states. A kinetic equation is solved analytically and numerically; the 
results agree well with the recent observations of Lewandowski et al. 
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In the last few years, the study of ultracold gases has 
generated a wealth of very interesting results. Spectac- 
ular examples are given by Bose condensed gases, but 
gases above their degeneracy temperature also provide 
exciting and unexpected results. For instance, recent ex- 
periments by Lewandowski et al. jjj have shown the exis- 
tence of a remarkable phenomenon, observed when a RF 
pulse is applied to a 87 Rb gas with two internal states, 
cooled by laser irradiation and evaporative cooling (but 
not Bose condensed). Since the two internal states are 
similar to two different species of atoms, the authors de- 
scribe their observation as a "segregation" between the 
species. They also mention that the differential Stern- 
Gerlach force, due to the magnetic gradient acting on the 
species, is too small to explain the segregation, which is 
actually related to interactions between the atoms. The 
purpose of the present article is to show that the "iden- 
tical spin rotation effect" (ISRE) provides a qualitative 
and quantitative explanation of the observations. 

The ISRE was introduced in || as a microscopic phe- 
nomenon taking place during a binary collision between 
two identical atoms with internal degrees of freedom, nu- 
clear spins for instance. The effect is a consequence of 
quantum indistinguishability; it introduces a rotation of 
each spin around their sum (in opposite directions for 
bosons and fermions) . For instance, if a single atom with 
a spin polarization in a given direction crosses a gas of 
identical atoms polarized in another direction, the spin of 
the transmitted atom undergoes a rotation; this is simi- 
lar to the rotation of the polarization of photons in the 
Faraday effect. On a macroscopic scale, the effect can 
affect transport properties of gases with internal states. 
For instance, ref. || considers a gas which is in a "classi- 
cal" regime in terms of equilibrium properties, but where 
quantum effects are important in binary collisions. When 
the density is sufficient to reach an hydrodynamic regime, 
this work shows the existence of transverse spin waves, 
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analogous to spin waves in degenerate liquid 3 He [|] . Sim- 
ilar predictions had been made independently by Bashkin 
]5| from a more macroscopic point of view, based on the 
notion of "molecular field" (or mean field) - see also the 
work of Levy and Ruckenstein || . Transverse spin waves 
in gases were subsequently observed in HJ, as well as 
in helium j|, ^] . 

Another prediction made in ref. || (end of section 1) is 
that the ISRE can also create "longitudinal oscillations" 
when a ir/2 pulse is applied to the sample, provided the 
transverse spin polarization is inhomogeneous. Here we 
show that the phenomenon described in ref. Q] is pre- 
cisely this effect, transposed to the pseudo spin associated 
with the two hyperfine levels relevant in the experiment, 
as foreseen by the authors who mention a "longitudinal 
spin effect" in their conclusion. The major differences arc 
that the experiment was performed at a density where the 
hydrodynamic regime is not reached, and that the spin 
oscillations are not of small amplitude. 

A point which emerged from the early studies on spin 
waves in gases, sometimes after vivid controversy, is that 
the effect of binary collisions in a gas are well described 
by a simple mean field calculation, provided one consid- 
ers forward scattering only. In condensed matter, each 
particle interacts at the same time with several others; it 
seems natural that their individual effects should be well 
averaged by the test particle so that mean field theory 
should apply. By contrast, in a dilute gas, particles are 
"free almost all the time" ; they interact only during brief 
collisions, with a single partner with which they can de- 
velop strong correlations. Indeed, in atomic physics, one 
rarely studies collision processes within mean field the- 
ory! Nevertheless, it turns out that the average effect of 
many collisions in the forward direction is equivalent to 
the results of mean field theory, if one replaces the real 
binary interaction potential by a pseudopotential involv- 
ing directly the scattering length (using the real potential 
would lead to meaningless results); the equivalence holds 
in the limit of low collision energy (the ISRE in the for- 
ward direction dominates over lateral scattering at very 
low energies since the corresponding "cross section" r^ d ~ 
diverges Q). 

In the experimental conditions of ref. the atoms 
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are in an axially symmetric magnetic trap elongated in 
the Ox direction. Initially the gas is at equilibrium with 
only state 1 populated. One then applies a 7r/2 RF 
pulse which, suddenly, puts all the atoms into the same 
coherent superposition of states 1 and 2, corresponding 
to a uniform transverse polarization of the pseudo spin. 
The system is then left free to evolve, and one observes 
the time evolution of the local densities n\ and n 2 . 

We begin with a qualitative physical discussion of the 
sequence of events. Since the field gradient creates an 
inhomogeneous spin precession, the gas develops a gra- 
dient of transverse spin orientation: correlations are cre- 
ated between position and transverse spin orientation. 
The free thermal motion of the atoms then creates cor- 
relations between velocity and transverse spin. Thus, a 
particle moving with a given velocity at point x gets a 
spin polarization which is not parallel to the average lo- 
cal spin polarization, so that the ISRE precession takes 
place. This makes its spin polarization leave the trans- 
verse plane and develop a non-zero value of its longitudi- 
nal component, with an opposite sign for different signs 
of the x-component of the velocity of the atom. The ap- 
pearance of this component indicates the beginning of 
an internal conversion, which eventually results in spa- 
tial separation of the atoms in different internal states. 
We emphasize that the apparent segregation is not the 
result of a spatial separation of atoms in fixed internal 
states, as for two different chemical species; on the con- 
trary, without changing their spatial position, the ISRE 
transfers atoms into internal states that depend on their 
motion. 

For a more quantitative discussion, we use a trans- 
port equation in terms of a time t dependent operator 
jo(r, p,t), which depends on position r and momentum 
p; pis the Wigner transform with respect to orbital vari- 
ables of the single particle density matrix; it remains a 
2x2 operator in the space of internal variables, corre- 
sponding to states 1 and 2. Instead of using the 4 matrix 
elements of p, it is often convenient to replace them by 
a local density / and (pseudo) spin density M in phase 
space defined by: 
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/(r,p,i)J + M(r,p,i) -a 
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where / is the unit operator in spin space and a the spin 
operator whose three components are the Pauli matrices. 
The kinetic equation for p(r, p, t) (see for instance H) is: 



dtp 
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where the second term is the usual drift term (777 is the 
mass of the particles); the third term (anticommutator) 
is the force term including both the effect of the trap- 
ping potential and of the mean field created by the other 
atoms; the fourth term (commutator) is a spin precession 
term containing the ISRE as well as some other contri- 
butions that we discuss below - this commutator is the 



term on which we focus our attention in this article. In 
the right hand side, the collision integral I co u [p] describes 
"real" collisions (lateral scattering as opposed to forward 
scattering, already included in the mean field); it can be 
obtained for instance from the LL transport equation [0 , 
or even take a more detailed expression containing "non- 
local collision terms" with r and p gradients, as discussed 
e.g. in the appendix of jlO) and pH . In fact, if we are 
mostly interested in a Knudsen regime, the precise ex- 
pression of I co u[p] is not needed. The effective potential 
U(r,t) is the spin operator: 



U(r,t) = U Q (r,t)I + U(r,t)-a, 
where the scalar component is defined by: 
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(7o-^-^+ g 22 n 2 +9lin 1 + 9 -f{n 2 +n 1 ). (4) 

Here V\ and V 2 are the external trapping potentials act- 
ing on states 1 and 2; the g e 's have the following ex- 
pressions in terms of the usual "coupling constants" g, 
proportional to the appropriate scattering lengths asso- 
ciated with the various possibilities for pair interactions 
between atoms in levels 1 or 2 [n2| : 



#11,22 = .911,22(1 + e)/2 ; g\ 2 = g d + eg t , 



(5) 



where g^ and g t refer to the direct and transfer process 
for two atoms in different levels. The number densities of 
atoms in levels 1 and 2 are 771,2; e = +1 (—1) for bosons 
(fermions). The vectorial component of U(r,t) is: 

U(r,t) = ^Me ||+ ^m(r,t); (6) 

e|| is the unit vector in the longitudinal spin direction, 
and 0(r, t) is: 

Ml = V 2 - Vi + 2 5 | 2 n 2 - 2^711 + 2g e 12 ( ni - n 2 ), (7) 
where the total density 77 and spin polarization m are: 

n(r, i) = J d 3 p f(r, p, t) ; m(r, t) = J d 3 p M(r, p, t) 

(8) 

and 77^2 = (n T m ||)/2- The first contribution (f7j) to 
U acts as a "local magnetic field" ; its average value over 
the sample can be removed in a uniformly rotating frame. 
The second contribution originates from the ISRE, and 
is proportional to the local spin polarization m (only m 
enters the ISRE commutator because the 1/k divergence 
of rl™ d ' at low fc's compensates the relative velocity 
factor of the collision integral). The commutator makes 
M precess around the momentum integrated local spin 
polarization; it does not affect the evolution of m itself, 
but can change the evolution of M for each value of p. 

A few simplifying assumptions are appropriate in the 
experimental conditions of [fij. The confining energy is 
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of order fc^T ~ 13 kHzx/i which is much larger than 
the mean- field interaction energy gn(0) ~ 140 Hzx/i and 
the differential trapping energy V\ — V% ~ 10 Hzxh. In 
the anti-commutator this allows us to retain only the 
confining energy of the harmonic trap: 

u * = \m [«V + coUy 2 + z 2 )] ; 0) 

here u> and uj ra( j are the axial and radial trapping frequen- 
cies. In the commutator of (||), Uq disappears, and only 
the vectorial component U plays a role. The cigar shaped 
trap has an axial frequency w/27r = 7 Hz, much smaller 
than the radial frequency uj ra d/2n = 230 Hz, so that the 
system is quasi one dimensional along the x axis. Assum- 
ing that radial local equilibrium is quickly established in 
the yz plane, we introduce on-axis value p(x,p,t) (inte- 
grated over radial momenta) . When averaged over radial 
coordinates and momenta, equation (o) becomes 



U(x, t) = Kl(x)e\\ /2 + eg{ 2 m(x, t) /4 



(10) 



where = V 2 — Vi + (g 22 — 3ii) n /2; note that the 
coupling constants are renormalized by a factor 1/2 upon 
averaging Q; we have assumed that 2g\ 2 — g{i + 322- 

With these assumptions the initial equilibrium 
Maxwcll-Boltzmann distribution f(x,p) solves the ki- 
netic equation (|), so that the dynamics after the pulse 
can be expressed in terms of M only. When lengths 
are measured in units of xt — \/kBT/muj 2 , momenta 
in units of px = ^rnksT and times in units of l/u>, 
equation (^|) simplifies into: 

2U M - M eq 

d t M+pd x M-xd p M- — x M~ (11) 

nuj lotm 

where a simple relaxation-time approximation has been 
made for I co u with a single parameter tm of the order of 
the time between collisions; the local equilibrium value 
M eq = m(x,£)exp(-p 2 /2)/\/27r. This is the equation 
that we now discuss. 

For small times we can use a time expansion 
M(x,p,t) = M(°! + M^f + M^ 2 h 2 /2 + .. and solve 
( |l l| ) to each time order. The spin distribution imme- 
diately after the tt/2 pulse is unchanged, except that 
M(°)(x,p) is perpendicular to ey. The density profile 
n(x) remains Gaussian, so that the Bohr frequency f2(x) 
does not vary in time. The effect of "real" collisions 
(right-hand side of (|Tl])) is neglected since we are inter- 
ested in small time behavior only. The result of this 
calculation is that m^(x,t) starts as t 4 with 
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eg{ 2 n(x) O (x)n(x) + f2 (x){n'(x) — 111(1)} /2H in di- 

mensionless form. Using the fact that the density n(x) 
remains Gaussian and restoring the units, we get: 



m||(x, t) n 2 — n\ 
n{x) n 



eg{ 2 n f2 x\ — 2Q x (ujt) 4 




time (ms) 
200 30CT 



400 
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FIG. 1: Time evolution of the spin polarization m at the cen- 
ter of the trap; my corresponds to the population difference 
between the two states; m_i_,l and mx,2 are the two compo- 
nents of the transverse spin polarization. The center-to-edge 
difference in Bohr frequency 8Q./2n is 12 Hz. 



The first factor in the right hand side is of the order of 
the dimensionless ISRE constant gi 2 n(0) / hui . Since its 
value in the experiment of ref. 13] is ~ 20, it is not surpris- 
ing that species separation could be observed in a time 
smaller than the trap period. The second factor involves 
the spatial variations of f2(x); near the center of the trap, 
a positive curvature implies a positive my (712 > ni), in 
accordance with the results of jj] (r/^ d ' < 0). Inserting 
the values of the parameters of the experiment of [Q into 
( |l2| ) leads to significant species separation (mil ~ n) for 
~ 25 ms, to be compared with the observed 30 — 50 ms . 

The maximum of the phenomenon can also readily be 
understood. For short times, we have seen that positive 
velocities along x correspond to one sign for the trans- 
verse orientation, and conversely. For times greater than 
2tt / VujSQ (5fl is the variation of f2(x) between the center 
and the edge of the cloud), both velocity signs become 
correlated to all spin directions in the transverse plane, 
so that the apparent segregation effect averages out to 
zero. Typical values taken from the experiment of ]lj] 
give ~ 100 ms for the maximum of the phenomenon. 

When, eventually, two separated species recombine un- 
der the effect of the restoring force of the trap, the ISRE 
plays no role anymore. The reason is merely that the op- 
erator associated with transverse spin is diagonal in the 
position representation (but not in the spin space), so 
that it can have non-zero value only if the wave-packets 
associated with each internal state overlap. In the ab- 
sence of transverse polarization, the system is equivalent 
to a classical mixture of two gases. 

Numerically, equation ( pd| ) can be solved by propagat- 
ing the initial distribution in time with the Lax-Wendroff 
method (see e.g. with parameters taken from ref. 

M. The Bohr frequency Q(x) is taken to be an inverted 
Gaussian of depth SQ and half- width xt- The dimension- 
less ISRE constant is gi 2 n(0)/huj — 20. The relaxation 
time is tm ~ 10 ms Q], so that ujtm — 0.3. The time 
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(a) (b) (c) (d) 



FIG. 2: Evolution of the particle density in state 1; 8SI/2-K 
and g 12 n(0)/hbj are (a) -12 Hz, 20 (b) Hz, 20 (c) 12 Hz, 20 
(d) 30 Hz, 30. Note the good agreement with Fig. 3 of Q, 
including the sign of the effect and "higher order effects" . 

evolution of the spin polarization at the center of the trap 
m(0, t) is shown in Fig. |[ with no adjustable parameter. 
The longitudinal spin polarization rises as predicted by 
(|l2|), reaches a maximum around 90 ms, and then oscil- 
lates and decays to almost zero after 300 ms. Even in the 
pure Knudsen regime (t~m = oo), a strong maximum of 
to 1 1 is reached around 100 ms. Fig. [ljalso shows how the 
other components of the polarization oscillate and decay. 
An interesting feature of the experiment of ref. Q is that 
neither the hydrodynamic nor the collisionless regime are 
valid along the axis as lotm ~ 1, so that a study of the 
full phase-space dynamics is necessary. 

Fig. [2] shows the time evolution of rii as a function of 
5f2, the variation of Q(x) between the center and the edge 
of the atomic cloud, and of gi 2 n(0) / ha; . When the cur- 
vature is zero (column (b)), no state separation occurs. 
Columns (a) and (c) show the effect of spin conversion 
for negative and positive curvature Stt (at the center). 



For negative curvature, the atoms in state 1 are pulled 
towards the center of the trap, whereas for positive cur- 
vature they are expelled from it. Column (d) exhibits 
what authors of ref. [|] call "higher order effects" , for 
sufficiently large values of <50 and gi2n(0)/tkJ. These fig- 
ures are in good qualitative agreement with Fig. 3 in ref. 

In conclusion, the ISRE plays an important role in the 
dynamics of cold gases with internal states. In the exper- 
iment of ref. |jj , this effect creates large longitudinal spin 
oscillations in a non-hydrodynamic regime. Our calcula- 
tions are also valid for fermions jl5) , where similar effects 
could be observed, in a case where g{ ± = g\ 2 = and 
the ISRE changes sign. Another interesting possibility 
is tuning the effect by changing g\ 2 at a Feshbach reso- 
nance [p.6| . Finally we note that, strictly speaking, our 
study is limited to non-degenerate gases; nevertheless, 
for non-condensed systems, most of the effect of degen- 
eracy can be included by simply replacing the Maxwell- 
Boltzmann distribution by the appropriate quantum dis- 
tribution p7| , so that no dramatic change is expected. 

Note: while this article was being written, we became 
aware of the work of Oktel and Levitov |l^] , who reach 
conclusions similar to ours. Nevertheless, they use an 
hydrodynamic expression for the evolution of the spin 
current, while here we put more emphasis on the inter- 
mediate and Knudsen regimes. 
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